Phase behavior of a nematic liquid crystal in contact with a chemically and 
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A nematic liquid crystal in contact with a grating surface possessing an alternating stripe pattern 
of locally homeotropic and planar anchoring is studied within the Frank-Oseen model. The combi- 
nation of both chemical and geometrical surface pattern leads to rich phase diagrams, involving a 
homeotropic, a planar, and a tilted nematic texture. The effect of the groove depth and the anchor- 
ing strengths on the location and the order of phase transitions between different nematic textures 
is studied. A zenithally bistable nematic device is investigated by confining a nematic liquid crystal 
between the patterned grating surface and a flat substrate with strong homeotropic anchoring. 
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PACS numbers: 64.70.Md, 61.30.Cz 



I. INTRODUCTION 



The contact of nematic liquid crystals (NLCs) with 
structured solid substrates offers the possibility to ma- 
nipulate the orientation of the adjacent liquid crystal 
molecules in a controlled way. For example, substrates 
with topographic structures influence anchoring of liquid 
crystal molecules and induce elastic distortions and flcx- 
oelectric polarizations within a contacting NLC. Align- 
ment of liquid crystal molecules can also be governed by 
the pattern of boundary lines between different regions 
on a flat substrate and the NLC elasticity. Adjacent re- 
gions of planar and homeotropic anchoring within a sin- 
gle substrate have been created, for example, by micro- 
contact printed polar and apolar thiols, respectively, on 
an obliquely evaporated ultra-thin gold layer. Whereas 
many experimental and theoretical studies have focussed 
on the understanding of the interactions and phase be- 
havior of NLCs in contact with ezi/ier geometrically struc- 
tured substrates (see e.g., [HEHHSESlEHIIl 
[lIEI El) or chemically patterned substrates (see e.g ., 

Hi ill in in m m, m m m mm in in mi), 

NLCs near geometrically structured and chemically pat- 
terned substrates have not been investigated yet. In view 
of the rich behavior of NLCs even on geometrically struc- 
tured or chemically patterned substrates, a combination 
of both surface treatments may open new possibilities for 
an improved performance of NLC cells. Here we study a 
NLC in contact with a grating surface possessing an al- 
ternating stripe pattern of locally homeotropic and pla- 
nar anchoring within the Frank-Oseen model [3fl [3l| . 
whereby the anchoring energy function is given by the 
Rapini-Papoular expression 32]. Taking into account 
both chemical and geometrical surface pattern is partic- 
ularly interesting because of the possibility of bistable 
anchoring, which is important for low power consump- 
tion liquid crystal displays. In bistable nematic devices 
there are two stable nematic director orientations which 



have substantially different tilt angles. The zenithally 
bistable nematic devices that have been studied recently 
consist of a NLC confined between a grating surface and 
a flat sur face El El El El- These surfaces are coated 
with a homeotropic agent. While such a device is char- 
acterized by both a high shock stability and a low power 
consumption, the deep grating structure leads to a re- 
duction of the contrast of the display EH ■ O n the basis 
of our calculations we expect that an additional chemical 
surface pattern allows one to use rather a shallow grating 
surface instead of a deep one. 

The structure of the paper is as follows. In Section ILTI 
two models of a patterned grating surface are introduced 
and phase diagrams of a NLC in contact with the model 
surfaces are discussed. In particular we show how phase 
transitions between different nematic textures vary with 
the groove depth and the anchoring strengths. The phase 
behavior and director profile of a NLC confined between 
a patterned grating surface and a flat surface subject to 
strong homeotropic anchoring is investigated in Sec. II I II 
Our results are summarized in Sec. II VI 



II. NLC IN CONTACT WITH A CHEMICALLY 
AND GEOMETRICALLY PATTERNED 
SUBSTRATE 

We first consider a semi-infinite system consisting of 
a NLC in contact with a single, patterned grating sur- 
face. Two models of the surface grating are investigated. 
In model A, a sinusoidal grating surface is assumed. 
This leads to the Euler-Lagrange equation for a two- 
dimensional director field, which is then solved numeri- 
cally. In model B, we choose a special form of the surface 
grating which, albeit not given explicitly by a simple for- 
mula, has the great advantage that the problem can be 
reformulated in terms of a semi-infinite system with a flat 
patterned substrate, and for the latter we have developed 
previously an efficient method of finding the equilibrium 
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[30l l3l| within the one-elastic-constant approximation, 
and the second term is the surface free energy (F s ) 
adopting the Rapini-Papoular form [32J. The anchor- 
ing strength is specified by a periodic step function: 
w(x) = ~u>h and w(x) — wp for values of x on the 
homeotropic and planar anchoring stripes, respectively 
(see Fig. ^| . Equation completely specifies the free 
energy functional for the system under consideration. 
Minimization of F[9] with respect to 0(x,z) leads to the 
Laplace equation with the following boundary condition 
at z = zq(x): 



FIG. 1: The system under consideration consists of a nematic 
liquid crystal (NLC) in contact with a grating surface with 
an alternating stripe pattern of locally homeotropic anchor- 
ing (white bars) and homogeneous planar anchoring (black 
bars). The stripes are orientated along the y axis perpendic- 
ular to the plane of the figure. The projection of the widths 
of the stripes onto the x-axis and the anchoring strengths are 
designated as pp = p/4, pp = p/4 and Wh, Wp, respectively, 
and the midlines of the planar anchoring stripes are located 
at x = np/2 with n = 0, ±1, ±2, ■ ■ ■ ; p is the period of the 
chemical pattern which is half the period of the surface grat- 
ing p. The tilt of the nematic director with respect to the 
z-axis is denoted by 8(x, z). 



director field 21). 



A. Model A: Sinusoidal grating surface 

Let us consider a NLC in contact with a single grating 
surface and assume that the surface profile zq(x) is given 
by 



Zq(x) = A sin (qx) 



(1) 



where A is the groove depth and p = 2tt /q is the period. 
Such a surface structure can be fabricated, for instance, 
using an interferometric exposure technique (see Ref. |l2^ 
and references therein). As Figure illustrates, the sur- 
face exhibits a pattern consisting of alternating stripes of 
locally homeotropic and homogeneous planar anchoring. 
The system is translationally invariant in the y direction. 
We focus on shallow grooves for which the nematic direc- 
tor n can be constrained to the plane perpendicular to 
the direction of the surface grooves. The orientation of 
n is therefore given solely by the polar angle 6(x, z) [see 
Fig. n . The free energy functional of the NLC reads 

F[9] = F d [9] + F s [6 ] 

K f p f°° 
= — dx dz [V9(x, z)] 2 

2 Jo Jz (x) 

1 r dx w{x) [- sm(9 (x))z' (x) + cos(fl (x))] 2 (2) 



v/i + («)) 2 

where 9q(x) = 9(x, zq(x)). The first term on the right 
hand side of Eq. (J2J is the distortion free energy (Fd) 



2K 



lim [-z' (x)d x + d z ] 9(x, z) 

-^z (x) 



w(x) 



2z' {x) cos(26» (x))} 



{[\-{z' a {x)f] sin(20 o (a;)) 



(3) 



and the second boundary condition: lirn^oo d z 9{x, z) = 
0. We solve this equation numerically on a sufficiently 
fine two-dimensional (x, z) grid. 

The periodic surface induces a certain anchoring di- 
rection (i.e., the orientation of n far from the surface) 



(e//) 



to 



which we call the effective anchoring direction 9 a 
distinguish it from the local anchoring directions of dif- 
ferent regions forming the surface pattern. In the ab- 
sence of external fields or competing surfaces n adopts 
the orientation 9 a e ^^ when the distance from the sur- 
face is large compared to the periodicity of the surface 
pattern. As we know .21], in the case of a flat sur- 
face (z'q(x) = 0) uniform nematic textures can exist: 
homeotropic (9(x, z) = 0) and planar (9(x,z) = vr/2), 
whereas it is clear from Eq. Q that no uniform texture 
(9(x,z) = const) is possible if z' (x) ^ 0. Nevertheless, 



the homeotropic (H) texture, defined by 9„ 



(eff) 



0, 



planar (P) texture, defined by 9a 



(eff) 



tt/2, can still 



exist even though 6{x, z) ^ const close to the surface. 
In other words, the presence of surface grating does not 
necessarily imply a tilted (T) nematic texture, corre- 
sponding to < 9 a e ^^ < tt/2. The existence of solu- 
tions corresponding to the H and P textures for arbi- 
trary groove depth follows from the assumed symmetry 
of the surface profile and anchoring strength function: 
zq(x) — zq(p/2 — x) and w(x) = w(p/2 — x). Since 
Eq. © is invariant with respect to the transformation 
x — > p/2 — x if 9{p/2 — x, z) = nir — 9(x, z), where n — 
or n = 1, we have 9 a e ^^ — for the antisymmetric so- 
lution (n = 0) and 9 a eff) = tt/2, for n = 1. Note that 
the above argumentation holds only if the functions Zq(x) 
and w(x) exhibit the same symmetry, although even in 
this case the tilted texture can be more stable than the H 
texture or P texture. If zq(x) and w(x) are not in phase 
only the tilted texture is possible. 
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B. Model B: Special form of surface grating 

Here we consider the surface profile zo (x) defined in an 
implicit form by 

zq(x) — Asm(qx) exp(— qzo(x)) . (4) 

It can also be expressed explicitly as 

M x ) = « _1 W (Aqsm(qx)) , (5) 

where Wq(C) is the principal branch of the Lambert W- 
function W(£) [35|. For small values of Aq, z (x) ob- 
tained from Eq. JSJ is very close to the sinusoidal grating 
given by Eq. QJ. An example of z (x) for A/p — 0.03 is 
shown in Fig. together with the sinusoidal grating sur- 
face used in model A. For both models, the minima and 
maxima of Zq(x) occur at x m i n — (n — 0.25)p and x max — 
(n + 0.25)p, respectively, where n = 0, ±1, ±2, • • • , but 
for model B \z (x m in)\ > \zo(x max )\- Since W(C) has a 
second-order branch point at £ = — e" 1 corresponding to 
Wo = — 1 36] our considerations are limited to groove 
depths smaller than A c /p = (27re) _1 w 0.06; for A — > A c 
we have z' (x m i n ) — > oc. 

Then we apply the following conformal mapping: 

t = x + Acos(qx) exp(— qz) , (6) 
u = z — A sm(qx) exp(— qz) (7) 

to express 9 as a function of (t,u). The condition 
z = zq(x) corresponds to u = in Eq. (0, thus, the 
distortion free energy expressed in terms of the new vari- 
ables has exactly the same form as for the flat substrate 
(cf. Eq. ©), i.e., 

K f p f°° 

Fd[0] = — dt du [(d t 9(t,u)) 2 + (d u 9(t,u)) 2 ] . 
1 Jo Jo 

(8) 

For the surface free energy we have 

1 f w(x)[-sm(0 o (t))z o (x)+co S (9 o (t))] 2 

*.m- 2 J Q (l-^ oW )(l + (4(^)2)3/2 > 

where 9o(t) = Q(t, u = 0), and x — x(t) is obtained from 
Eqs. © and for u = 0. 

By comparing Eq. J5J) with Eqs. (JSJ and © we see that 
due to the conformal mapping the semi-infinite system 
with the rough patterned surface has been mapped onto 
the semi-infinite system with a flat patterned surface and 
a modified form of the surface free energy. The Green's 
function method can now be used to express 6(t, u) that 
satisfies the Laplace equation in terms of the boundary 
function 9o(t), which allows to treat F as a functional 
of 9o H3- Thus, instead of solving the Laplace equation 
with suitable boundary conditions (see Eq. 10) in the 
two-dimensional xz plane we minimize F[9o] numerically 
on a one-dimensional grid of variable t. 
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FIG. 2: Comparison of the surface profiles given by Eq. JTJ 
[solid line] and Eq. @ [dashed line] for A/p = 0.03. 



C. Phase diagram 

Figure [3] displays the phase diagram constructed as 
a function of the groove depth A and the homcotropic 
anchoring strength wh, for three values of the planar 
anchoring strength wp. The textures P, T and H are 
separated by the lines of phase transitions which we re- 
fer to as the P-T, T-H and P-H transition, respectively. 
Both model A and B predict that the P-T transition is 
always continuous, whereas the T-H transition can be ei- 
ther first- or second-order, depending on the ratio A/p. 
We note, however, that the direct first-order P-H transi- 
tion, which occurs when the grooves are sufficiently deep, 
follows only from model A. 

In the limit of small wh the planar texture is sta- 
ble. When wh increases, and the grooves are shallow, 
the second-order P-T transition occurs. Our calcula- 
tions show that for a fixed value of wp the location of 
this transition is rather independent of the groove depth, 
and when wp increases it moves towards larger values of 
wh- For deeper grooves, the line of the continuous P-T 
transition terminates in a critical end point on the first- 
order T-H/P-H transition line. In the case of flat surface 
(A = 0), the T-H transition is always continuous and can 
exist only for pwp/K < 9. If pwp/K > 9 the H texture 
does not exist even for large values of wh- For grating 
surfaces, the T-H transition changes from the second- to 
first-order when the groove depth increases. The change 
occurs at the tricritical point whose position shifts to 
higher values of A upon decreasing wp. 

We emphasize that in our model of the NLC in contact 
with a chemically patterned but flat substrate (A = 0) 
the phase transitions between the P or H texture and 
the T texture are always continuous, and there is no di- 
rect transition between the H and P textures. On the 
other hand, we observe a first-order P-H transition for a 
chemically uniform sinusoidal surface but for rather deep 
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FIG. 3: Phase diagram of a NLC in contact with a single sinu- 
soidal grating surface [model A, see Eq. Q and Fig. with 
alternating stripe pattern of locally planar and homeotropic 
anchoring as a function of groove depth A and the strength of 
homeotropic anchoring wh for three values of the strength of 
planar anchoring wp. The solid and the dashed lines denote 
first- and second-order transitions, respectively, between pla- 
nar (P), tilted (T), and homeotropic (H) nematic textures. 
The solid and open circles mark tricritical and critical end 
points, respectively. For comparison, the results obtained 
from model B [see Eq. (QJ] are shown with the same line code, 
but for A < 0.057 p. 

grooves (A > 0.25 p), similarly to the results obtained for 
asymmetric surface grating structures 0- However, as it 
is apparent from Fig. |3 a combination of both chemical 
and geometrical surface pattern can lead to the first-order 
H-T transition, and also reduce the value of A/p above 
which the H-P transition can be found. The repercus- 
sions of this result on possible bistable nematic devices 
are discussed in Sec. III. Note also that since the nematic 
texture induced by a periodic substrate is characterized 
by one of the effective anchoring directions: B a e ' = 0, 
6»i e//) = tt/2 or < 9 { a ff} < tt/2, the observed phase 
transitions can be regarded as anchoring transitions. 

In Fig.[3]we have also compared the predictions of mod- 
els A and B. We observe a good agreement for small val- 
ues of A/p, whereas for deeper grooves some deviations 



FIG. 4: Phase diagram in the (pwn/K, pwp/K) plane for: 
(a) A < A , (b) A = A ~ 0.075 p and (c) A > A . The 
meaning of symbols and lines is the same as in Fig. 01 Dia- 
gram (a) follows from the calculations based on model B for 
A/p — 0.02. In the case of flat (A = 0) chemically patterned 
substrate (not shown here) both phase transitions are contin- 
uous everywhere (the tricritical point on the H-T line escapes 
to oo). Diagrams (b) and (c) correspond to model A but they 
are only schematic. At A — Ao the tricritical point and the 
critical end point replace each other. 



appear. This reflects the fact that the difference between 
the grating profiles defined by Eqs. and J3J becomes 
more pronounced with increasing groove depth. Since 
the application of model B is limited to rather small val- 
ues of A/p, for the reasons already discussed, this model 
is unable to describe the first-order transition between 
the homeotropic and planar textures. However, the ad- 
vantage of model B is that the precise determination of 
phase diagrams for small ratios A/p becomes feasible, al- 
lowing one to study various aspects of the phase diagrams 
in more detail and with considerably less computational 
effort than within model A. Thus, to obtain the over- 
all picture of the phase transitions studied we have used 
information from both models. 

In Fig. we present a sequence of phase diagrams 
in the (pwn/K, pwp/K) plane, for three values of the 
groove depth: (a) A < A , (b) A — A and (c) A > A , 
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FIG. 5: Three-dimensional phase diagram in the 
(pwp/K, pwh/K, A/p) space. The meaning of sym- 
bols and lines is the same as in Fig. |3] The solid and dashed 
lines correspond to intersections with planes A = const, and 
the grey area corresponds to the surface of the first order 
transitions. Note the evolution of the tricritical point and 
the critical end point when A/p changes. 



where Aq /p w 0.075. In all diagrams the P-T transition is 
continuous. In case (a), there is a tricritical point on the 
H-T transition line. For A — > its position tends to the 
infinite value of wh (see Fig. [5} , which is consistent with 
the flat substrate case. Upon increasing A the tricriti- 
cal point moves towards (pwn/K = 0, pwp/K — 0) and 
eventually, at A = Aq (case (b)), it reaches the origin, 
and the H-T transition becomes first-order everywhere. 
In case (c), the continuous P-T transition line does not 
extend to the origin, as in cases (a) and (b), but termi- 
nates at the first order transition line, which for A > Aq 
consists of two pieces: the H-T and P-H lines. Note that 
when A passes Aq the tricritical point is replaced by the 
critical end point or vice versa. 

To better visualize the relation between phase dia- 
grams presented in Figs. |21 and 0] a schematic three- 
dimensional phase diagram in the space spanned by 
(pwp/K, pwn/K, A/p) is also shown (see Fig. (SJ. The 
two-dimensional plots can be deduced from Fig. [S] by 
making intersections parallel to the plane (A/p, pwn / K) 
or (pwu/K, pwp/K). 

Finally we note that one can also consider the situation 
when homeotropic stripes are on the sides and planar 
stripes are on the hills (and in the pits) of the surface 
profile (see Fig.QJ. In such a situation, the H-T transition 
line can be easily deduced from the P-T transition line of 
the above-discussed case by interchanging wh with wp 



FIG. 6: Phase diagram of a NLC confined between a pat- 
terned sinusoidal grating surface (see Fig.0 and a flat surface 
with strong homeotropic anchoring as a function of groove 
depth A and the strength of homeotropic anchoring wh for 
two values of the strength of planar anchoring mp. The solid 
and dashed lines denote first- and second-order transitions, re- 
spectively, between a homeotropic (H) and hybrid aligned ne- 
matic (HAN) texture. The solid circles mark tricritical points. 
The width of the cell is D/p = 2. 

and ph with pp; the same holds for the line of the P-T 
phase transition. 



III. ZENITHALLY BISTABLE NEMATIC 
DEVICE 

We now turn our attention to a NLC confined between 
a chemically patterned sinusoidal surface (which corre- 
sponds to model .4.) and a flat substrate with strong 
homeotropic anchoring 6(x, D) = 0, where D is the thick- 
ness of the cell. In the following calculations we keep 
D fixed, but vary both the anchoring strengths and the 
groove depth at the grating surface. The numerical so- 
lution of the Euler-Lagrange equation demonstrates the 
existence of two nematic director configurations, namely 
the homeotropic (H) texture, where the director field is 
almost uniform and perpendicular to the flat surface, and 
the hybrid aligned nematic (HAN) texture, where the di- 
rector field varies from homeotropic to nearly planar ori- 
entation through the cell (c.f. Fig.[7|bclow). Of course no 
planar texture, which exists in the semi-infinite case, can 
be found because of the strong homeotropic anchoring on 
the flat substrate. 

Figure displays the phase diagram plotted as a func- 
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FIG. 7: The tilt of the nematic director 8(x = 0, z) with re- 
spect to the z-axis for a NLC confined between the patterned 
sinusoidal grating surface shown in Fig. and a flat surface 
with strong homeotropic anchoring. The dashed and the solid 
lines represent the profile for the coexisting H and HAN tex- 
tures for the planar anchoring strength pwp/K = 2.45 and 
groove depth A/p — 0.075; see also the schematic nematic di- 
rector profiles on the right side bar. The flat surface is located 
at z = D = 2p while the grating surface is at zo(x — 0) = 
(see Fig. [Q. The location of the hills of the grating surface is 
marked by the dotted line. 



tion of the groove depth A and the homeotropic anchor- 
ing strength wr on the grating substrate for two values 
of the planar anchoring strength wp and for a fixed cell 
thickness D. For weak homeotropic anchoring the HAN 
texture is stable provided the groove depth is smaller 
than A up p. For deeper grooves with A > A upp no HAN 
texture is found because distortions of the nematic direc- 
tor field are too costly in the presence of the dominat- 
ing homeotropic anchoring. Upon increasing the groove 
depth, the HAN-H transition changes from second- to 
first-order at the tricritical point. 

The tilt angle 9(x, z) in the middle of the planar stripe 
(x = 0) is shown in Fig. [7] as a function of z for the coex- 
isting H and HAN textures. The nematic director is al- 
most parallel to the z-axis in the H texture (for z > O.Qp), 
while it changes smoothly from 9(0, z = D) = to 
8(0, z = 0) < 7r/2 in the HAN texture. In the latter case 
9(0, z) decreases roughly linearly as the distance from 
the flat surface increases. Moreover, the numerical cal- 
culations exhibit that the nematic director field is rather 
independent of x except of a thin surface layer (of width 
« p/(2n)) at the grating surface. 

Finally, we emphasize that the groove depths for which 
zenithal bistability is found within the present model are 
considerably smaller than those considered theoretically 
(A = 0.8 p) [HEl and used experimentally (A = 0.35 p) 
[l2j for a NLC confined between a pu re g eometrically pat- 
terned surface and a flat substrate 0, 0] . The calcu- 
lated value of the average tilt of the nematic director just 
above the grating surface in the HAN texture is compara- 
ble to the one obtained in the case of the deeper grooves 

mm 



IV. SUMMARY 

We have applied the Frank-Oseen model together with 
the Rapini-Papoular surface free energy to a NLC in 
contact with a single grating surface possessing an al- 
ternating stripe pattern of locally homeotropic and pla- 
nar anchoring (Fig. ^| or confined between the patterned 
grating surface and a flat substrate. Phase diagrams and 
nematic director profiles are determined numerically with 
the following main results. 

(1) A homeotropic, a planar, and a tilted nematic tex- 
ture have been found for the NLC in contact with a single 
patterned grating surface (Figs. 0-0- Both second- and 
first-order transitions between the tilted texture and the 
homeotropic or planar texture are possible. Furthermore, 
for appropriate values of the groove depth and local an- 
choring strengths one can also observe a first-order tran- 
sition between the homeotropic and planar textures. It is 
worthwhile to emphasize that the combination of chem- 
ical and geometrical surface pattern reduces the groove 
depth above which this transition occurs as compared to 
a pure geometrically patterned surface. 

(2) For the NLC confined between a chemically pat- 
terned sinusoidal grating and a flat substrate which in- 
duces strong homeotropic anchoring, we have determined 
a first-order phase transition between a homeotropic tex- 
ture and a hybrid aligned nematic texture for rather 
shallow grooves (Fig. 0). In the homeotropic texture, 
the nematic director is almost uniform and perpendicu- 
lar to the flat surface, while the director field varies from 
homeotropic to nearly planar orientation in the hybrid 
aligned nematic texture (Fig.[7J. Building on the results 
shown in Figs. [S] and [7| it seems possible to achieve a 
zcnithally bistable nematic device without the use of a 
deep surface grating. 

It is worthwhile to emphasize that not any combination 
of the chemical pattern and surface grating is capable of 
inducing first order transitions between the homeotropic 
and tilted textures or the homeotropic and planar tex- 
tures in a semi-infinite NLC system. For instance, if the 
periods of the chemical pattern p and of the surface grat- 
ing p are the same then the surface roughness merely 
changes the location of the second-order transitions ob- 
served in the flat substrate case, at least in the range of 
the groove depth studied in this work. Moreover, our nu- 
merical calculations suggest that the optimum ratio p/p 
for the occurrence of the first order transitions in the case 
of relatively shallow grooves is p/p = 2 shown in Fig. ^ 

Finally, it is instructive to consider the situation when 
the chemical surface pattern is not exactly in phase with 
the surface grating profile. In such a case, for the reasons 
mentioned in Sec. Ill Al there are no true homeotropic 
and planar textures, hence the accompanied second-order 
phase transitions cease to exist. However, first order 
transitions between a pseudo-homeotropic texture and 
the tilted or hybrid aligned texture should still exist, 
where the term pseudo-homeotropic means that the tex- 
ture is characterized by w 0. Note also that the 
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tricritical point in the phase diagrams shown in Figs. [3]- 
[SJbecomes a critical point at which the difference between 
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